Dynamical behavior of generic quintessence potentials: constraints on key dark energy 

observables 
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We perform a comprehensive study of a class of dark energy models - scalar field models where 
the effective potential can be described by a polynomial series - exploring their dynamical behavior 
using the method of flow equations that has previously been applied to inflationary models. Using 
supernova, baryon oscillation, CMB and Hubble constant data, and an implicit theoretical prior 
imposed by the scalar field dynamics, we find that the ACDM model provides an excellent fit to 
the data. Constraints on the generic scalar field potential parameters are presented, along with the 
reconstructed w(z) histories consistent with the data and the theoretical prior. We propose and 
pursue computationally feasible algorithms to obtain estimates of the principal components of the 
equation of state, as well as parameters wo and w a . Further, we use the Monte Carlo Markov Chain 
machinery to simulate future data based on the Joint Dark Energy Mission, Planck and baryon 
acoustic oscillation surveys and find that the inverse area figure of merit improves nearly by an 
order of magnitude. Therefore, most scalar field models that are currently consistent with data 
can be potentially ruled out by future experiments. We also comment on the classification of dark 
energy models into "thawing"' and "freezing" in light of the more diverse evolution histories allowed 
by this general class of potentials. 



I. INTRODUCTION 

Is the universe accelerating because it is dominated by 
vacuum energy? The answer to this question has become 
a holy grail of cosmology; in addition to the increasing 
improvements in existing constraints [1-5] , there are mul- 
tiple ambitious experimental efforts currently planned to 
help find the answer. Vacuum energy would be spatially 
homogeneous and time- invariant, possessing a dark en- 
ergy equation of state, w = pde/pde, equal to —1 iden- 
tically and at all times. Finding robust evidence for a 
deviation from this prediction would be tremendously 
important and would suggest yet another cosmological 
mystery: it would indicate that dark energy is more com- 
plicated than the simplest model, Einstein's cosmological 
constant. 

Mapping out the history of the equation of state of DE 
(or dark energy density), and in particular any variation 
in redshift of w(z), is therefore one of the fundamen- 
tal goals in cosmology. The simplest way to do this is 
to measure a single parameter that describes the time- 
dependence of w(z), such as its derivative (e.g. [6-9]), 
while the most general approach is to directly reconstruct 
w(z) from the distance-redshift or expansion rate-redshift 
data [10-12]. With any given test, a difficult question 
arises: if we measure that the equation of state at some 
pivot redshift is consistent with —1 with some error, is 
it worth spending a large amount of time and resources 
in measuring the temporal variation of w — which, one 
could naively guess, would then be small or zero in any 
realistic model of dark energy? 

In this paper we address the following question: given 
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the current constraints on the equation of state of 
dark energy that are increasingly tightening around the 
ACDM value of —1, and assuming a well-defined class of 
models that allows diverse behavior in the dark energy 
sector, is it worth pursuing better and more expensive 
cosmological probes in hope of detecting any deviation 
from the ACDM value? Our analysis differs from a large 
body of previous work in that we approach the ques- 
tion in a very general way. Instead of considering very 
specific examples of DE models (say, specific power law 
scalar field potentials), we encompass all models within 
a prescribed framework, constrain this entire class with 
the currently available data, and address how future im- 
provements upon these constraints will affect our ability 
to test the ACDM paradigm for dark energy. To perform 
these tasks, we use the formalism of dynamical trajecto- 
ries in model parameter space that was previously applied 
to slow-roll inflation models. 

We make no attempt to study the individual or com- 
parative merits of specific cosmological probes; nor do we 
consider optimal survey design for any particular probe. 
Such studies have been carried out by various authors in 
the past (e.g., [13-22]), and many others have obtained 
constraints on standard descriptions of the dark energy 
sector [23-34] as well as various parametric descriptions 
of the energy density [35-46]. Moreover, comparison of 
cosmological probes strongly depends on the correct char- 
acterization of the systematic errors. Instead, we concen- 
trate on asking specific questions about a large, precisely 
defined physical class of dark energy models using a com- 
pilation of current and (simulated) future data. In the 
process we propose several recipes for converting between 
different relevant parameter bases, and also for describ- 
ing the cosmological data. In this sense, we significantly 
extend the work of Sahlen et al. [47] who adopted a simi- 
lar approach of Monte Carlo reconstruction of the scalar 
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field dark energy models. 

The paper is organized as follows. In Sec. II we de- 
scribe the motivation for this work and briefly describe 
the class of dark energy models that we consider. In 
Sec. Ill we lay out the equations necessary to compute 
the dark energy history of each model, and outline our 
main assumptions that further define our framework. In 
Sec. IV we describe the Markov Chain formalism we 
adopt, specify the initial conditions, and discuss extrap- 
olation between the low redshift and the CMB era. In 
Sec. V we define some of the derived parameters that de- 
scribe the dark energy history, such as the principal com- 
ponents and wo, w a and w p i vo t. In Sec. VI we present 
the cosmological constraints on the various parameters 
and functions, while in Sec. VII we discuss the cosmo- 
logical implications and figures of merit. We conclude in 
Sec. VIII. Appendices A and B describe in detail, respec- 
tively, the current and future cosmological data we have 
used, as well as the likelihoods assigned to them. 



II. MOTIVATION AND FRAMEWORK 

We would like to "scan" through a variety of dark en- 
ergy models with a time-dependent equation of state in 
order to explore the generic allowed ranges of dark en- 
ergy sector parameters at redshifts where the observa- 
tions are powerful, and then test how well ACDM can be 
distinguished from dark energy models in terms of cur- 
rent and future constraints. In most general sense, scan- 
ning through dark energy models is impossible as a single 
physical framework for dark energy is currently nonexis- 
tent. Therefore, we have to specialize to one of the several 
classes of dark energy models, or else describe the back- 
ground evolution of the universe without recourse to a 
physical model of dark energy. Here we adopt the former 
approach and choose perhaps the most widely considered 
model: a rolling scalar field, or quintessence [48-54]. To 
scan through the space of scalar field models, we adapt 
the Monte Carlo reconstruction formalism that has pre- 
viously been applied to inflationary dynamics. Note that 
the application of the flow equations to the DE case re- 
quires some modifications compared to the case of in- 
flation. Even though the acceleration of the universe 
implies that we recently entered an inflationary period, 
the densities of other components such as matter are not 
negligible, and this prevents us from relating the Hubble 
parameter to the slow-roll parameters in an easy way as 
with inflation. In particular, so-called Hubble slow roll 
parameters (e.g., [55, 56]) cannot be used since the Hub- 
ble parameter is not changing very slowly because of the 
matter component. In fact, there are no small parame- 
ters in the DE flow equation formalism, as the potential 
slow-roll parameters are not necessarily small. Even for 
a steep potential with large "slow-roll" parameters much 
greater than unity, Hubble friction may still enable slow 
roll with small kinetic energy of the field and w(z) ~ — 1. 
Therefore, the allowed ranges for the analogous param- 



eters in the dark energy case are determined entirely by 
the data. 



III. METHODOLOGY 

A. Scalar held equations 

We start with the Klein-Gordon equation for the single 
scalar field 



' + 3^+^=0, 



(1) 



where the overdot represents a time-derivative. We im- 
pose two requirements on the rolling field: 

1. The field is not allowed to turn around during its 
roll on the potential (i.e. models where d<j> changes 
sign are rejected). If the turnaround happens, the 
model is assigned zero likelihood. 

2. At initial time, the field is only allowed to roll down 
the potential slope. 

These requirements are in general not required; how- 
ever they simplify the initial assumptions while not weak- 
ening the generality of the framework. Fields with 
rapidly changing sign of dej) — violating (1) — describe 
potentials with the time-averaged equation of state near 
zero (e.g., an axionic field), therefore not describing the 
DE. Dark energy models therefore obey the first require- 
ment. The requirement (2) was imposed to speed up 
scanning of the models. We have explicitly considered 
fields initially rolling up the potential, and found that 
nearly all of them turn around, thereby violating require- 
ment (1). 

We can rewrite Eq. (1) in terms of dark energy's equa- 
tion of state and energy density as a function of scale 
factor 



( 2 ) 

where <f> = 4>/m p \. Hereafter, we consider all dimension- 
ed quantities in units of m p i and drop the tilde for sim- 
plicity. 

We integrate Eq. (2) starting at some redshift z start . 
To do so, we apply initial conditions as values of the 
equation of state and energy density of the field at z s t a rt , 
w start = w(z sta . Tt ) and Og| rt = fbE^start)- Using the 
Friedman equation and the definitions of w and Ode, we 
find the expressions for w and f^E at any given time as 
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and where a start is scale factor at the initial time and V a = 
V(<p — 0) is the potential at initial time of integration. 
We need one further initial condition: 
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Integrating Eq. (2) with initial conditions w start ^ Slp| rt 
and Eq. (7), together with relations (3-6), gives the full 
dark energy history out to arbitrary redshift. In this 
work, we integrate the equation starting at z stai t = 3. 
While we can obtain constraints on the parameters of in- 
terest at any given redshift, we choose to consider them 
at z = which is close to the redshift where the sensi- 
tivity of observations is greatest. We have checked that 
the constraints at z = do not change appreciably if we 
start integrating either at z s tart = 2 or z s t ar t = 5. 

Eq. (2) can be integrated forward or backward in time. 
We have explicitly checked that the forward and back- 
ward integration give precisely the same dark energy his- 
tory (and values of all parameters as a function of red- 
shift) provided that the final conditions of integration 
forward in time were used as initial conditions of the in- 
tegration backward in time, or vice versa. While is it 
not explicitly needed to solve the system of equations 
specified above, we will also need the following quantity 
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B. The flow equations as a potential generator 

The one extra ingredient we require in order to incor- 
porate this formalism into cosmological parameter esti- 
mation is a potential generator which parameterizes the 
potential V((f>) in as general a manner as possible [57]. In 
order to accomplish this task we adapt the inflationary 
flow equation formalism. Since it was first proposed a few 
years ago [56, 58], this formalism has been used to gener- 
ate a large number of inflationary models in a relatively 
model-independent way for the purpose of Monte-Carlo 
reconstruction (e.g. [59, 60]). Its principal advantage 
is that it does not rely on any specific particle physics 
model, but rather starts with the slow-roll formalism and 
provides an approximate expansion of the effective poten- 
tial in a hierarchy of slow-roll parameters. It has recently 
been incorporated directly into cosmological parameter 



estimation in a way minimizes the effects of the lack of 
knowledge of the measure of initial conditions [61, 62], 
and we follow the spirit of this latter work. 

The dark energy case differs from the inflationary case 
because it is not described by slow roll parameters. How- 
ever we can still specify an infinite hierarchy of parame- 
ters in terms of the derivatives of the potential: 



l X = 




t > 1, 



(9) 



where prime denotes derivatives with respect to <j>, and 
later we will refer to : A as rj and 2 A as £. These parame- 
ters are related to each other through an infinite system 
of coupled first order differential equations. The trajec- 
tory specified by this infinite system is exact; however, in 
practice we must truncate it at some finite order, and ob- 
tain an approximate potential. Truncating the hierarchy 
of flow parameters at the term M A means that M+1 A = 
at all times as well. From Eq. (9), it also follows that 
d (M+2) v y d( ^(M+2) = o at all times. This truncated hi- 
erarchy is therefore closed and has an analytic solution 
[63] ; it simply describes a polynomial of order M + 1 in 
V(4>) with 



V(4>) = V [1 + A 1( f> + A 2 4> 2 + ... + A 
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The coefficients A,, with i > 1, are written in terms of 
the starting values of the flow parameters as 
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where A\ — — \/l67re start specifies the direction the field 
is rolling initially (an opposite choice of this direction 
and A\ will lead to identical results provided the initial 
velocity in Eq. (7) is given a minus sign). The expansion 
is in field space, not redshift space, and without loss of 
generality we assign (f> = at the starting redshift, z sta rt- 
As you move away from 2 s tart, different values corre- 
spond to different redshifts for different models; thus the 
expansion is not coeval between models. Further, the 
specification of the parameters f A star t at z star t specifies 
the potential over all redshifts. Constraints obtained can 
be translated to any desired redshift, and therefore the 
formalism does not contain an additional assumption of 
a pivot-redshift. 

In the inflationary context, the fact that the flow pa- 
rameters are small (because the inflaton has to roll slowly 
over ~ 60 e-folds of inflation) means that one can make 
an argument for retaining terms to a high order in the 
hierarchy (say M + 1 = 10), assuming that successive 
higher order flow parameters can be chosen from nar- 
rower and narrower priors. That is, a slow roll expansion 
like Eq. (10) must describe any generic smooth, flat po- 
tential quite accurately over a wide range of the field roll. 
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TABLE I: The classes of dark energy models are shown in the 
first column. The second and third columns list the numbers 
and names of the parameters describing the effective poten- 
tial. Note that proxies for the initial field value and initial 
potential energy of the field are parameters not listed in this 
table — ng| r * and w start , respectively. 



Model 


M+l 


Potential Parameters 


2 V{4>) parameters 


2 


,start ^start 
t i 7 


3 V{4>) parameters 


3 


start start tstart 



In the dark energy formalism, on the other hand, the 
flow parameters are not necessarily small any more, and 
one can no longer make the assumption that successive 
higher order flow parameters can be chosen from nar- 
rower and narrower priors - i.e. the potential is not nec- 
essarily described by a slow roll expansion. However the 
dark energy field does not need to roll for a large number 
of e-folds in order to explain the observed acceleration 
(unlike in inflation), and therefore Eq. (10) only needs 
to describe the field over a narrow range in its evolution. 
Also, there are no currently conceived cosmological obser- 
vations that can provide us with information to constrain 
more than a handful of potential parameters. Therefore 
we don't necessarily need to make a slow roll expansion 
valid over a wide range of field evolution. 

Instead, Eq. (10), truncated at a low order, now serves 
as a generator of potentials, where the values of the flow 
parameters at £ s tart are generated in a Monte-Carlo fash- 
ion to simulate large numbers of potentials valid over the 
redshift range [0, z sta rt]- We explicitly check below, us- 
ing wide ranges for the priors on the initial values of 
the potential parameters (i.e. making no assumptions of 
smallness), that the w(z) histories generated by truncat- 
ing the hierarchy at M+l = 2, M + 1 = 3 and M+l = 5 
are generically qualitatively similar, and therefore our pa- 
rameter estimation is not very sensitive to the order of 
this truncation. We posit that this similarity is due to 
two reasons: firstly because the w(z) history only cap- 
tures a very limited amount of information about the de- 
tailed shape of the potential, and secondly because such 
dysmorphic potentials have a greater chance of the scalar 
field changing direction, thus causing those models to be 
excluded due to our theoretical prior. 

In the parameter estimation we carry out below, we 
choose to close the hierarchy at M + 1 = 2 and M+l = 3 
and consider the models defined in Table I. This is equiv- 
alent to setting (</>) and all higher derivatives to zero 
in the first case, and setting (</>) and all higher deriva- 
tives to zero in the second. Finally, note that proxies for 
the initial potential energy of the field (equivalently, pa- 
rameter Vb) and initial field velocity <fo are the parame- 
ters Og| rt and w start j respectively. 



IV. MARKOV CHAIN MONTE CARLO 
ANALYSIS 

We use a Markov Chain Monte Carlo (MCMC) tech- 
nique [64-69] to evaluate the likelihood function of model 
parameters. The MCMC is used to simulate observations 
from the posterior distribution V(a\x), of a set of param- 
eters given event x, obtained via Bayes' Theorem, 



V{0\x) = 



V(x\6)V(6) 

jv(x\o)v(e)de' 



(12) 



where V{x\6) is the likelihood of event x given the model 
parameters and V(0) is the prior probability density. 
The MCMC generates random draws (i.e. simulations) 
from the posterior distribution that are a "fair" sample 
of the likelihood surface. From this sample, we can esti- 
mate all of the quantities of interest about the posterior 
distribution (mean, variance, confidence levels). A prop- 
erly derived and implemented MCMC draws from the 
joint posterior density V(6\x) once it has converged to 
the stationary distribution. 

We use four chains per model and a conservative 
Gelman-Rubin convergence criterion [70] to determine 
when the chains have converged to the stationary dis- 
tribution. For our application, 6 denotes a set of cos- 
mological parameters: the starting value of the energy 
density and equation of state of dark energy, f2f)§; rt and 
ui start , and the starting values of the potential parame- 
ters as specified in Table I. This set of parameters fully 
determines the dark energy history from z — z star t un- 
til today. In addition, we marginalize over the angular 
size of the acoustic horizon 9a (a proxy for the Hubble 
constant H a ) and the fractional physical energy density 
in baryons, ttbh 2 as nuisance parameters. The universe 
is assumed to be to be flat; thus the fractional matter 
density is given by fi m = 1 — Q\. 

The full details of the cosmological data used, along 
with their likelihood functions, are given in Appendix 
A. Essentially, we use a combination of supernova data 
(Supernova Legacy Survey [5]), baryon oscillation re- 
sults from the Sloan Digital Sky Survey [71], cosmic mi- 
crowave background constraints from the WMAP exper- 
iment [26], and the Hubble Key Project [72] measure- 
ment of the Hubble constant. Appendix B details the 
projected future cosmological data, where we assume the 
same cosmological probes but with the expected, smaller 
statistical errors (as well as estimates of the systematic 
errors). 

Note that we do not know a priori the ranges of the ini- 
tial conditions, since experimental measurements probe 
a specific temporal average of ^de and w. Therefore, we 
try to keep the initial conditions as general as possible, 
and fully marginalize over the maximally uninformative 
flat priors, 
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FIG. 1: Left panels: representative histories w(z) for a number of randomly chosen sets of initial conditions, shown for 
illustration without applying cosmological data constraints. Models which are freezing (see Sec. VII B) are color-coded in blue, 
while the models color-coded in black are neither thawing nor freezing. We find that thawing models are very rare. Right 
panels: the corresponding phase plots in the w — dw/dlna plane, showing the evolution of these parameters as a function of 
redshift. The evolution histories originate at z s tart = 3. The orange diamonds represent z — 0. The parameterization assumes 
2-parameter (top row), 3-parameter (middle row) and 5-parameter (bottom row) descriptions of V(<f>). 
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e start G [0,oo] 

^start g [.^oq] 

^tart g [.qo^o]. ( 13 ) 

Note that we had to impose an upper limit on £ start in 
order to get the chains to converge: the upper limit in 
^start j g essentially unconstrained, both for current and 
future data, and our imposed limit assures convergence 
without affecting the results. We then apply the following 
Monte-Carlo algorithm: 

1. Pick random initial values from the priors in 
Eq. (13), together with values of 9a and fi^/i 2 . 

2. Integrate the Klein-Gordon equation; Eq. (2) to 
get f^DE at the present epoch and w(z). Compute 
the derived parameters fl m h 2 and Hq, and also the 
principal components and wq and w a (see next Sec- 
tion). 

3. Compute the likelihood of this model (using the 
likelihood functions defined in the Appendices). 

4. Move to the new model using the Metropolis- 
Hastings sampler. 

5. Repeat from step 2. 

Figure 1 shows randomly chosen DE histories w(z) as- 
suming 2, 3 and 5-parameter descriptions of V(<j>) (left 
panels; top, middle, and bottom rows respectively). For 
illustration, we show here the results of Monte Carlo 
draws from the prior, before applying cosmological con- 
straints. To a first approximation, the three cases look 
qualitatively very similar. In particular, both show ex- 
amples of freezing models — where w{z) is asymptoti- 
cally approaching —1 (blue curves), using the language of 
Caldwell and Linder [73], while we see very few thawing 
models — where w(z) is asymptotically receding from 
— 1. We also see examples of models that are neither 
purely freezing nor thawing (black curves); we comment 
on this in more detail in Sec. VII B. The right panels of 
Figure 1 show the corresponding phase plots (i.e. the evo- 
lution of each of the models above on the w — dw/d In a 
plane). These plots therefore show typical histories of 
w(z) which are generated by the prior on the parameters 
which are going to be subsequently confronted with the 
data using the MCMC method. 

Note that we have cases where w(z sta rt) > 0. While 
these models may in principle not be ruled out by the 
low-z data (SNe, BAO), they may well not be viable if 
the energy density in DE is high enough at z > z stai t 
to spoil the successful formation of structure. A strong 
test of such models is imposed by the CMB measure- 
ment of the distance to recombination, encapsulated by 
the angular size of the acoustic horizon 9 a- To apply 
the CMB acoustic scale test, it is clear that we need to 
know w(z) at z > z s tart5 however this is difficult as it 
does not seem reasonable to extrapolate our local ex- 
pansion of V(4>) to the early universe and integrate the 



equation for 4> backward to z w 1000. Instead, we choose 
a simpler extrapolation in the equation of state, and as- 
sume that the equation of state is constant beyond our 
starting redshift, w(z > z start ) = w(z stai t)- While some- 
what ad hoc, this choice seems reasonable, as the only 
way that most models with w(z > z star t) > could sur- 
vive the data cut is if their potential had a feature, so 
that w{z) at high z becomes again sufficiently negative. 
We implicitly ignore this small class of potentials that is 
unconstrained by the data. Further, we have repeated 
our analysis without the 9a constraint (that is, without 
the constraint from distance to the last scattering sur- 
face) and found small differences at z = z star t — dark 
energy models with w(z sta rt) > are allowed without 
the 9a essentially because the early history of dark en- 
ergy is completely unconstrained in that case. However, 
the low-redshift constraints are nearly unchanged with 
or without 9a imposed. Models that would be ruled out 
with the 9a (that is, those with the "incorrect" distance 
to the last scattering surface) are often ruled out by the 
SNc-t-BAO+CMB+i/o combination alone because they 
do not revert to being close to ACDM at low redshifts 
where the constraining power of the data is the greatest. 
We conclude that our ansatz w(z > z s tart) = w(z sta rt) 
has no discernible effects on the final low-redshift con- 
straints. 

The following sections detail the dark energy observ- 
ables on which we obtain constraints. 

V. DERIVED PARAMETERS: DEFINITIONS 

In what follows we adopt the 2-parameter description 
of V(4>) (see Eq. (10)) as our fiducial case, and also con- 
sider the 3-parameter description as a test of the sensi- 
tivity of our results to the truncation of V (</>). Therefore, 
we have a total of at least four parameters that describe 
the background cosmology: at least two describing the 
potential ( e start and ry start ) ; and then Q^e* and w start . 
Note that we will be mainly interested in constraints on 
the present-day values of the first four parameters, which 
we denote with the superscript "0" . Finally, we will not 
be very interested in the constraints on the remaining 
parameters, 9a, Hq, Q m h 2 and flbh 2 , because we found 
that, as expected, the constraints on these parameters 
largely reproduce the prior given to them by the obser- 
vational measurements (the details of which are in the 
Appendices). 

In addition to presenting constraints on our fiducial 
parameters and looking at dark energy histories via w(z), 
we would like to report errors on the commonly used 2- 
parameter description of the equation of state (wo,w a ) 
where [7, 74] 

w(z)=w +w a — — . (14) 

1 + 2! 

However, as seen in Figure 1, the scalar field models 
evidently do not follow the relation (14) over the redshift 
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range considered, though they (presumably) do over the 
range of redshifts strongly probed by the data. Therefore, 
we need an algorithm to assign the best-fitting wq and w a 
to a given dark energy history. We do that in two steps: 
we first determine the best-fitting principal components 
computed from an idealized Fisher matrix analysis, then 
convert the first two of them into the familiar equation 
of state parameters. We now describe this procedure in 
detail. 



A. Principal components 

A convenient way to describe a given dark energy 
model is to compute the best-measured principal compo- 
nents (PCs) of w(z) for each dark energy model history 
[75]. This is a form of data compression, where many pa- 
rameters needed to describe the dark energy sector are 
replaced by a few that describe the quantities that are 
best measured by a given set of cosmological probes. The 
main advantage of the PCs is that they are model in- 
dependent — in other words, they are independent of 
parameterizations of dark energy density or equation of 
state. The weight of the best-measured principal compo- 
nent is a model-independent predictor of what redshift 
range is best probed by a given cosmological probe - 
the first PC peaks at z ~ 0.2 for future SNe measure- 
ments [75, 76], z ~ 0.5 for weak lensing [77, 78] and 
z ~ 0.7 for baryon oscillations. The shapes of principal 
components also depend, albeit more weakly, on the sur- 
vey specifications and the fiducial cosmological model. 

In order to compute the principal components from a 
given survey, one needs to compute joint constraints on a 
large number of parameters that determine the function 
w(z), then diagonalize their covariance matrix. While 
relatively straightforward to do in the Fisher matrix for- 
malism and with simulated data [75-77, 79-81], this task 
is very challenging with actual data as accurate computa- 
tion of the covariance matrix of a large number (~ 10-50) 
of highly correlated parameters is required, although the 
computational requirements are significantly less severe 
if one only wants a rough resolution of the principal com- 
ponents in redshift as done in Refs. [46] and [82]. 

Note that we could in principle compute the exact prin- 
cipal components of our scalar field model equations of 
state, since we have the full w(z) of each model in our 
Markov chains; it is just a matter of outputting these 
quantities at a number of redshifts, computing their co- 
variance matrix from the information in the chains, and 
diagonalizing it. However, the equation of state values 
at different redshifts in scalar field models are highly cor- 
related since the scalar field paradigm together with the 
few-parameter description of the effective potential mod- 
els only allows specific dark energy histories. These cor- 
relations render the principal components of the scalar 
field equation of state very different from the usual PCs 
of a general w(z). While it may be of interest to work 
with the scalar field PCs, this is largely terra incognita, 
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FIG. 2: The first three principal components that are best 
measured by the current data (solid, dashed and dotted 
dark/blue line) and the future data (solid, dashed and dot- 
ted dark/blue line light/orange line). The components were 
computed starting with 40 piecewise constant values of the 
equation of state linearly distributed in the redshift range 
< z < 2 and using the Fisher matrix formalism, as de- 
scribed in [75]. Note that the future PCs have peak weight 
at somewhat higher redshifts, which directly reflects the ex- 
tended redshift reach of future SNe and BAO observations. 



and in this work we instead concentrate on the usual PCs 
for a completely general w(z) computed for a combina- 
tion of surveys that we consider. We leave the study of 
scalar-field-specific PCs for future work. 

Therefore, we adopt a hybrid approach here: we com- 
pute the principal components assuming a Fisher matrix 
as an approximation to the inverse of the true covariance 
matrix, and also assuming theoretically unconstrained 
dark energy histories. We use the actual cosmological 
data, current or future (from Appendices A and B), ap- 
proximating the CMB information by the distance to the 
last scattering surface which has been marginalized over 
the sound horizon parameters; this analysis closely par- 
allels that in [75, 77]. We parameterize the equation of 
state in 40 piecewise constant values uniformly distribu- 
tion in < z < z max - We evaluate the original Fisher ma- 
trix, invert it, eliminate the parameters not correspond- 
ing to w(z) (such as Hq and f2 m ), and diagonalize the 
resulting matrix in the equation of state parameters to 
obtain the principal components. 

Figure 2 shows the shapes &i{z) of the first three prin- 
cipal components of the current data (solid, dashed and 
dotted dark/blue line) and the future data (solid, dashed 
and dotted light/orange line). The curves have been 
smoothed in redshift with a cubic spline, and are nor- 
malized arbitrarily in this plot. Note that the weights of 
the future PCs peak at somewhat higher redshifts, which 
directly reflects the extended redshift reach of the future 
SNe and BAO observations. 

Given a dark energy history w(z), the amplitude of zth 
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principal component, on, is given by [75] 



f 

on = 

Jo 



w(z) ei(z) dz 



(15) 



where the principal components are normalized so that 
ei(z) dz — 1, and the error in on is determined by 
the cosmological datasets used. For each model from 
the Markov chains, wc compute the principal component 
coefficients according to Eq. (15). Given that we are now 
using the pre-computed principal components that are 
obtained for a general w(z), the parameters cti will not 
be uncorrelated in our scalar field dark energy analysis. 



B. The equation of state parameters 

We would also like to show our constraints in terms of 
more familiar 2-parameter description of dark energy w 
and w a . Of course, our goal is not to simply fit (wo, w a ) 
to the cosmological data as often done in the literature — 
instead, we seek to describe each scalar field model that 
we generate with some effective (wo,w a ). Mapping the 
true redshift-dependent equation of state of a model into 
a finite number of parameters is not trivial. The most ob- 
vious way to proceed would be to fit (wo,w a ), together 
with other cosmological parameters, to each dark energy 
model in our chains. Unfortunately, this procedure would 
be too time-consuming, as we would need to perform a 
multi-parameter fit to each one of millions of dark en- 
ergy models that we generate in our MCMC procedure. 
Clearly, a faster way to construct (wo,w a ) is needed. 

Here we propose and adopt an alternative, extremely 
simple algorithm: we simply convert the first two prin- 
cipal components (011,0:2) into (wo,w a ). This approach 
is justified because the first two PCs carry essentially all 
the necessary information about the effects of dark en- 
ergy dynamics on the expansion of the universe over the 
observable scales. From Eqs. (14) and (15) it follows that 
we can define 



w 



w a = 



0=201 - ai02 

01-02 
ai — ct2 

01-02 



where 



0i 



J Bi(z) - 



+ z 



■ dz. 



(16) 
(17) 

(18) 



Equations (16) and (17) arc now our definitions of the 
parameters wo and w a , given a dark energy history w(z) 
which determines the on. We have explicitly checked on 
individual examples that the two-parameter equation of 
state closely follows the true w(z) over the redshift range 
most effectively probed by the data. Note too that the 



constraint wq > — 1, which follows from w(z) > — 1, is not 
strictly obeyed by wq obtained in this way since Wo and 
w a are now essentially a fit to the dark energy equation 
of state history. We will return to the efficacy of this 
approximation later in Sec. VI D. 



C. Equation of state pivot 

The pivot value of the equation of state parameter is 
the value of w{z) at the specific redshift where the equa- 
tion of state is best constrained and has a "waist" [14] . In 
other words, we require that w p i vo t = wo+Aw a and w a be 
uncorrelated, where A is a coefficient to be determined. 
In practice, having obtained the (wo,w a ) corresponding 
to a model for each step in the Markov chains, we can 
compute the covariance matrix of these two parameters 
and diagonalize it to extract (w p i vo t,w a ). It is easy to 
show that A = — Cov(w , w a )/Cov(w a , w a ) and 



Cov(w , w a ) 

Wpivot = W - -, r W a 

Cov(w a , w a ) 



(19) 



where Cov(x, y) stands for elements of the covariance ma- 
trix element. Although we make no specific use of the 
pivot redshift, for completeness, it is given by 



Zpivot 



Cov(w , w a ) 



Cov(w , w a ) + Cov(w a , w a ) ' 



(20) 



Now we will present our constraints on these dark en- 
ergy observables. 



VI. BASIC RESULTS 
A. Fundamental parameters 

In Fig. 3 we show the constraints from the current data 
on e, n, f^DE, and w at redshift zero (left panel; super- 
scripts "0") and redshift z star t = 3 (right panel; super- 
scripts "3"). In Fig. 4 we show the equivalent constraints 
for future data, assuming a ACDM fiducial model. 

The z = 3 constraints (right panel of Fig. 3) arc not 
particularly tight, as the only direct probe of this red- 
shift comes from the CMB. Nevertheless, we see that the 
equation of state is limited to be t<; start < 0; this con- 
straint is directly due to the CMB location of the first 
peak 9a, as models with w(z) > at high redshift are 
increasingly dark energy dominated and produce an in- 
correct distance to the last scattering surface. Note that 
the best-fit ACDM cosmology has Og| rt « 0.05. 

In contrast, constraints at z = (left panel of Fig. 3) 
are excellent. In particular the equation of state w° and 
the first slow-roll parameter e° arc determined with great 
accuracy. Our choice to show constraints at z = is 
somewhat arbitrary, as we would expect similarly good 
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FIG. 3: Solid lines show the marginalized 2D-joint 68% and 95% probability contours (off-diagonal panels) and ID marginalized 
probability distribution (diagonal panels) for the two potential parameters e and rj, as well as Que, and w. The color coding 
in the off-diagonal panels shows the marginalized probability density in these 2D parameter spaces, ranging from red for the 
highest density to blue for the lowest. The left and right sets of panels show constraints at z = and 2 s tart = 3 respectively, 
using the current cosmological data as described in Appendix A. We do not show the parameters Ho, Q m h 2 and Qbh 2 which 
largely reproduce the applied measurement priors. Note that the constraints at z s tart = 3 are considerably weaker than at 
z — 0, since the constraining power of the data is concentrated at lower redshift. 



constraints at z < 0.5 in general. Rather than these 
quantities at particular time slices, it is of more inter- 
est to show the reconstructed histories of dark energy, as 
well as quantities that reflect the time-averaged behav- 
ior of dark energy (e.g. wq, w a , principal components) 
as we do below. Nevertheless, we pause to comment 
that these constraints are complementary to the usual 
approach where the dark energy sector is parameter- 
ized via a single, constant equation of state parameter 
w. While neither one of our fundamental parameters 
in Fig. 3 is constant in time, the constraints on them 
are tight because of the combination of the theoretical 
prior and cosmological data. Constraints from these two 
different approaches, model-independent w = const and 
scalar field time-dependent w(z), do not necessarily need 
to be in agreement — nevertheless, we do expect many 
salient features to be the same, such as the preference for 
^de ~ 0.75 and w° « — 1, and this is reflected in the left 
panel of Fig. 3. 

One interesting point to note is that the expected fu- 
ture constraints on the potential parameters e and r\ im- 
prove to a significantly lesser degree than those on Ode 
and w when compared with the constraints from the cur- 
rent data. This suggests that our knowledge of the lo- 
cal shape of the effective potential will not improve very 
much in the future compared to that of the velocity and 



the potential energy of the field. 

In our formalism we can directly compute the field ex- 
cursion from z = 3 to z = for each model in our chains 
by integrating Eq. (2), and hence obtain constraints on 
it: we find A0/m pl = 0.09±0.03. Thus the field evolution 
in the range probed best by the data is small. Therefore 
it is not surprising to find that the constraints on the 
shape of the potential are not strong, since such a nar- 
row range of the potential is probed. In addition, unlike 
in the case of inflation, there is no theoretical slow-roll 
prior to help constrain the shape of the potential over a 
larger (f> range than that probed by observations. 



B. Principal components 

Figure 5 shows the 2D-joint 68% and 95% CL con- 
straints in the (0:1,012) (left panel) and (02,03) (right 
panel) planes. Constraints using the current data, de- 
scribed in Appendix A, are shown in blue/dark regions, 
while constraints from the future data, described in Ap- 
pendix B, are in orange/light regions. As mentioned ear- 
lier, the future data are centered on the ACDM model 
(a, = —1 for all i) by fiat. We overlay the actual con- 
straints with ~ 1000 models randomly selected from our 
prior without being subject to the data constraints. Note 
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FIG. 4: Same as Fig. 3, but parameter constraints are for future data as described in Appendix B, denned at z = (left) and 
Zstart = 3 (right). It is interesting to note that constraints on the potential parameters e and r\ improve to a significantly lesser 
degree than those on Que and w when compared with the constraints from the current data. 



that the future constraints are much better than the cur- 
rent ones (this is further discussed in Sec. VII A). While 
the density of the models itself depends on the prior, it is 
clear that both the current and the future data exclude 
a significant fraction of them. Finally, the marginalized 
constraints for the individual principal component pa- 
rameters are shown in Table II. 

The constraints on the principal components exhibit 
features that are interesting to comment on. Principal 
components of general dark energy models are by defini- 
tion uncorrelated, and this is clearly not the case from 
Fig. 5; the reason is, as explained in Sec. VA, that we 
have used the general PCs and constrained them assum- 
ing scalar field models. The principal components are 
therefore correlated, but this does not change any of the 
scientific results in the paper. 

Moreover, one can note four well-defined edges in the 
left panel of Fig. 5 for either dataset. The condition 
a i > —1 is always strictly obeyed since the first principal 
component is a non-negative weight over the equation of 
state and w(z) > — 1 for scalar fields; see Fig. 2. For the 
same reason the condition > —1 is obeyed, though not 
strictly — the second PC has a both positive and nega- 
tive weight and, for a small fraction of our w(z) models, 
this leads to «2 < —1 (a few of the dotted examples in 
the left panel of Fig. 5, around the value a\ ~ 0.7, ex- 
hibit this behavior). The long diagonal edge seen in this 
panel is also due to the requirement w(z) > —1 and its 
precise slope and intercept are determined by the shapes 
of the first two principal components. Finally, the fuzzy 



edge with roughly ot\ < —0.85 (for the current data) is 
due to the cosmological data that effectively restrict the 
weighted average of the equation of state to be less than 
this value. This last limit is more stringent for the future 
data, and dotted models from the prior, which have not 
been subjected to the data, are obviously not restricted 
in this direction. 



C. Parameters wq, w a and u> p i V ot 

Constraints in the wo-w a plane are shown in the left 
panel of Fig. 6. As in Fig. 5, we have shown the current 
and future constraints (the latter assuming ACDM ) as 
well as ~ 1000 individual models from the prior. More- 
over, there are two sharp edges that all models obey due 
to the requirement w(z) > —1; this is similar as with 
the principal components shown in Fig. 5. Likewise, the 
third, long side of the triangle contours are due to the cos- 
mological data which, to a first approximation, impose an 
upper limit to an average of the equation of state. For the 
current data, this constraint is roughly 6wq + w a < —5. 

Following the simple procedure outlined in Sec. VC, 
we can obtain constraints on the pivot value of the equa- 
tion of state. Constraints in the (w p i vo t,w a ) plane are 
shown in the right panel of Fig. 6. These constraints 
are similar to those on (wo,w a ) except that w p i vo t and 
w a are uncorrelated by the definition of w p i vo t- As with 
the principal component cvi, the relatively sharp edge at 
Wpivot = — 1 is due to the fact that w(z) > —1 for scalar 
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2 V(0) parameters 2 V(0) parameters 




FIG. 5: Left panel: 68% and 95% CL constraints from current data (larger/blue contours) and future data (smaller/orange 
contours) on the principal components ai and Qf2- Right panel: same constraints, but on the principal components «2 and 
03. The future data have been centered on the ACDM model (on = — 1 for all i), while the current data reflect the actual 
constraints. The points show ~ 1000 DE models simulated randomly from the priors, that passed our basic requirements set 
forth at the beginning of Sec. Ill but have not been subjected to the data. Note that the points lying within the region allowed 
by the future data constraint have been overplotted by the constraint for clarity — models from the prior go all the way down 
to the ACDM solution. 




FIG. 6: Left panel: 68% and 95% CL constraints from current data (larger/blue contours) and future data (smaller/orange 
contours) on the equation of state parameters wo and w a . Right panel: same, but the parameters are now w p i vo t and w a , 
where the pivot w p i vot is the value of w(z) at the specific redshift where wo and w a are decorrelated. As in the previous two 
figures, the future data have been centered on the ACDM model (coordinates ( — 1,0) in both panels); the current data reflect 
the actual constraints, and the points show ~ 1000 DE models simulated randomly from the priors. 



field models. Future constraints are about an order of 
magnitude better than the current ones, as noted above. 

Finally, Table II shows the constraints on wq, w a , and 
Wpivot- Note that wq is a two-tailed distribution while 
Wpivot is a significantly skewed one-tail distribution with 
a hard prior u> P i VO t > — 1 ; where the probability density 
abruptly falls at the boundary. This is because while wq 
is allowed to go below —1 due to our method for assign- 
ing (wo, w a ) for a given model w(z), decorrelating these 
parameters recovers the prior on the original evolution 
histories, w(z) > — 1. Further, contrary to the expecta- 



tion from the usual Fisher matrix analyses, the strong 
non-Gaussianity of the probability distribution function 
means that, while the if p i vo t and w a are decorrelated, it 
is not true any more that we have the best constraint on 
the equation of state at z p i VOt . However, it is still useful to 
go to the uncorrelated parameter basis (w p i vo t,w a ) even 
without this latter benefit, as w P i VO t is a better approx- 
imation than wq to the real w(z) over the range probed 
by the data, in particular recovering the theoretical prior 
WpWot > — 1. 
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D. Equation of state reconstruction 

Finally, in Fig. 7 we show the 68% and 95% CL regions 
on the reconstructed equation of state w(z). The recon- 
struction is straightforward, as the values of w(z) at a 
number of redshifts are written out for each of our mod- 
els as derived parameters in the MCMC. The left panel 
shows the constraints from the current data, while the 
right panel shows reconstruction from the future data. 
In the left panel, the dashed lines show the 68% and 95% 
reconstruction limits using the (w ,w a ) parametrization. 

The constraints are very good even for the current data 
set. They impose a sharp upper limit on the equation 
of state at low redshift, and progressively get weaker at 
higher z. Moreover, w(z) < is required at any red- 
shift since otherwise, dark energy would be increasing in 
importance relative to matter at early times and would 
spoil the distance to the last scattering surface (as well as 
structure formation). Finally, note that the constraints 
using the (wo, w a ) parametrization are in excellent agree- 
ment with the true constraints out to z w 0.7 — there- 
fore, our definition of wq and w a , Eqs. (16) and (17), 
does an excellent job in reconstructing the dark energy 
history out to this redshift. In fact, we have checked 
that the rms difference between the true and (u; , un- 
reconstructed distance out to z = 3 is only about 0.8% 
(for models from the posterior distribution for the cur- 
rent data), though it degrades to 3% by z — 1089. As 
discussed in Sec. VII D, the fitting accuracy to z = 3 is 
significantly better (by about a factor of 4) than the ac- 
curacy to which the z = 3 distance is determined from 
the current data, and therefore there is no bias in us- 
ing the (wo,w a ) fit in this redshift range; however it is 
dangerous to extrapolate the fit out to higher redshift. 

Note that this Monte Carlo reconstruction is very dif- 
ferent in spirit from the direct non-parametric recon- 
struction of the equation of state from distance measure- 
ments [10, 11, 35, 83-86] — the former is done via mod- 
els, scalar field in this case, while the latter is completely 
independent of any models (although the non-parametric 
reconstruction still requires parametric fits to smooth the 
noisy distance-redshift curve observed by, say SNe la). 
Heuristically, the Monte Carlo reconstruction follows the 
dark energy histories generated by the assumed models, 
and the models correlate the equation of state values at 
different redshifts. This is the reason that we obtain sig- 
nificantly better and more reliable constraints than in 
the general non-parametric case. Conversely, the recon- 
struction presented here is clearly model-dependent, and 
necessarily less general than in the non-parametric case. 
Fig. 8 shows explicit examples of w(z) histories taken 
from the posterior. 

E. Three-parameter V(cf>) constraints 

There is no reason to limit the effective potential to 
be of second order: our parametrization naturally allows 



potentials of arbitrary order, and we are only limited by 
how many parameters can be constrained with current 
or future data. We illustrate the constraints obtained 
with the polynomial of 3rd order in Fig. 10. The top 
left panel shows the constraints on e°, 77°, f^ E , w° and 
£° (i.e., at z = 0), while the top right panel shows the 
same constraints at 2; sta rt = 3. The middle row shows 
the constraints in the (wo,w a ) plane (left panel) and the 
(wpivot, Wa) plane (right panel), while the bottom row 
show the reconstructed w(z) for the current data (left 
panel) and future data (right panel). 

The three-parameter V{<j>) constraints on the five fun- 
damental parameters are somewhat weaker than those 
for the two-parameter case (shown in Fig. 3). In particu- 
lar, the three-parameter constraints seem to get progres- 
sively worse as we go from e° to 77° to £°. Nevertheless, 
the limits are still interesting, especially those on 
and w°. However, constraints on the parameters wq, w a , 
Wpivot and the principal components aj are essentially 
unchanged from the two-parameter case. This means 
that, despite the more complicated potential that is now 
harder to constrain, the range of dark energy histories 
constrained by the data is largely unchanged. This could 
also be seen by comparing the range of models present 
in the posterior distribution; see the left three panels in 
Fig. 1. Thus, we conclude that our constraints are only 
weakly sensitive to the truncation order of the polyno- 
mial series. 



TABLE II: 68% C.L. constraints on the first three principal 
components, and on parameters wo, w a , and w p i vo t. The 
second column shows the constraints from the current data, 
while the third column shows errors from the future data, 
assuming the underlying ACDM cosmology (for which all pa- 
rameters are —1 except w a which is zero). All rows report 
the usual two-sided error bars, except w p i vo t where we report 
one-sided error bars because of the mathematical constraint 
Wpivot > — 1 where the probability density abruptly falls at 
the boundary. Recall that the constraints are significantly 
better than in the general non-parametric case because of the 
implicit theoretical prior of using scalar fields. 



Parameter 


Current data 


Future data 


cti 


-0.949i™i 


fiduciall^^ 


a 2 


-0.8991°;^ 


fidudallffi 


a 3 


-0.9011™! 


fiduciallHii 


w 


-0.9781^1 


fiduciali^^ 


w a 


0.1971S 


fiducial!™*? 


^pivot 


< -0.957 


< -0.987 
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FIG. 7: Shaded regions show the 68% and 95% CL constraints on the equation of state w(z) from the current data (left panel) 
and future data (right panel). The future data assume that the fiducial cosmology is precisely ACDM w(z) = — 1. In the 
left panel, the dashed lines show the 68% and 95% reconstruction limits using the (wo,w a ) parametrization. Note that the 
approximate (wo,w a ) reconstruction is in excellent agreement with the actual constraints at z ^ 0.7. 



2 V(0) parameters (examples from the posterior) 3 V(0) parameters (examples from the posterior) 
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z z 



FIG. 8: Follows the format of Fig. 1, except these figures show the evolution histories of a number of models which were 
accepted steps in the MCMC. Therefore, whereas Fig. 1 shows examples of models generated by the prior, this figure shows 
examples of models in the posterior. 



F. (Lack of) sensitivity to initial conditions 

In section III A we have briefly mentioned that the low 
redshift constraints are largely independent of the start- 
ing redshift of integration z s tart- We further illustrate 
this statement in Fig. 9 where the shaded regions show 
the reconstruction of the equation of state w(z) for the 
current data as in the left panel of Fig. 7. However now 
the dashed lines show the reconstruction assuming the 
starting redshift z sta rt = 5 rather than z s tart = 3 (the 
epochs 3 < z < 5 are not shown for clarity) . Clearly, the 
constraints at z < 0.7 are essentially insensitive to the 
starting redshift. 

Solid lines in Fig. 9 show the reconstruction of w(z) 
also with Zstart = 5, but now forcing the scalar field to 



start at rest (so that w start — — 1). Clearly the field 
rolls slowly in the beginning as it gains speed, and the 
equation of state is close to —1. Nevertheless, we find 
again that memory of this transient behavior is lost at low 
redshift (z < 0.7) where the constraints are unchanged 
relative to the fiducial case with unconstrained w start . 

These two exercises suggest that the low redshift con- 
straints are sensitive only to the details of the potential 
and not initial conditions. However, we have decided to 
leave a full study of this issue, including imposing physi- 
cally motivated constraints on the initial conditions that 
take into account the field evolution at z > z s tart, for 
future work. Here we have adopted a completely em- 
pirical approach and allowed maximally general initial 
conditions of the scalar field. 
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FIG. 9: Shaded regions show the 68% and 95% CL constraints 
on the equation of state w(z) from the current data, and 
are identical to those in the left panel of Fig. 7. Dashed 
lines show the reconstruction assuming the starting redshift 
Zstart = 5 rather than z s tart = 3 (the epochs 3 < z < 5 are not 
shown for clarity). Solid lines show the reconstruction also 
with z s t a rt = 5, but now forcing the scalar field to start at 
rest (so that u; start = —1). Clearly, the constraints at z < 0.7 
are insensitive to changing the initial conditions. 



VII. COSMOLOGICAL IMPLICATIONS 

We now present some implications of our constraints 
regarding the figures of merit, the classification of the 
dynamics of our models, and the preference of ACDM 
over the more complicated scalar field models. To begin 
with, it is clear from all results presented so far that the 
ACDM model is a very good fit to the current data. This 
result is in full accord with previous analyses [25-28, 87]. 



A. Current vs. future data 

One of the questions we set out to answer with our 
approach was whether it is worth pursuing future experi- 
ments, given that the current constraint are increasingly 
converging on the ACDM model. A view of Fig. 5 shows 
that the answer is clearly affirmative. For example, the 
inverse area of the 2-dimensional 68% CL contour in the 
wo-Wa plane, sometimes taken as the figure of merit for 
the power of cosmological constraints [14, 88], is about 
8 times better for the future constraints than the cur- 
rent ones; it is 12 times better if we compare the 95% 
CL contour areas. For the first two principal compo- 
nents, the FoM is a similar factor of ~ 10 better for the 
future dataset. Clearly, if dark energy is observation- 
ally distinguishable from the ACDM model, upcoming 
surveys will have a significantly greater power to reveal 
this. Therefore, we eagerly expect the next generation of 
experiments, led by Planck, Joint Dark Energy Mission 
(JDEM) and Large Survey Telescope (LST). 



B. Thaw or freeze? 

There has recently been a surge of interest in clas- 
sification of dark energy models. In particular, it has 
been emphasized that well known dark energy models — 
scalar field models with potentials that are power laws 
or exponential functions of the field — naturally fall in 
classes of "thawing" or "freezing" [73, 89-91] depend- 
ing on whether they are asymptotically receding from or 
approaching the state of zero kinetic energy where the 
equation of state is — l 1 . Note that freezing dark energy 
models have been studied at least as far back as [48], 
while thawing models date back at least to [50]. 

With the benefit of our framework we can quantita- 
tively answer questions about the division of models. Our 
parameterization is more general than that in most pre- 
vious studies as we consider all quadratic/cubic polyno- 
mials for V(cj>) with maximally uninformative priors for 
the initial speed and energy density of the field. On the 
other hand, we concentrate mainly on the low redshift 
universe and do not attempt to model the effective po- 
tential during the "dark ages" of the universe and earlier. 
Therefore it is interesting to compare our findings with 
those from previous studies. 

We define a given model as "thawing" if dw/dlna > 
uniformly, and "freezing" if dw j d In a < uniformly, and 
as neither thawing nor freezing if dw /din a changes sign, 
in the interval z s tart = 3 to z = 0. From our prior alone, 
with no data cut applied (as in left panels Fig. 1) we find 
that about 40% models are freezing, 0.4% are thawing, 
and about 60% are neither. If we consider the models in 
our chains, i.e. in the posterior distribution (as in Fig. 8), 
we find that about 74% are freezing, a negligible fraction 
(0.05%) are thawing, and 26% fall into neither category. 
Hence the thawing fraction, small to begin with, is nearly 
completely eliminated once the data have been applied. 

One might consider that, since thawing models have 
ui start ~ — 1 and SlpE * ^ 1 an< ^ we nave uniform pri- 
ors in those quantities, that the prior is biased against 
generating these models. However, it is not necessarily 
true that the posterior will therefore be biased against 
them. If the data clearly favored thawing models over 
freezing models, the MCMC would "learn" to concen- 
trate very close to the initial conditions that favor these 
models. Instead we see that the fraction of freezing mod- 
els in the posterior increases over that of the prior, and 
the already-small fraction of thawing models becomes yet 
smaller after the application of cosmological data, which 
favor models that have w(z) ~ —1 at low redshifts. 

This situation is analogous to the case of Monte Carlo 



1 Note that Ref. [73] excludes consideration of potentials with a 
flat direction, V'((p) = 0, on the grounds that these possess a hid- 
den cosmological constant; more general quintessence potentials 
allowing for flat directions have been considered in the literature, 
eg. [92]. 
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FIG. 10: Constraints obtained assuming V(cj>) is a polynomial of 3rd order. The top left panel shows the constraints on e°, rf, 
Qde, io° and £° (i.e., at z — 0), while the top right panel shows the same parameters evaluated at z s tart = 3. The middle row 
shows the constraints in the (wo,w a ) plane (left panel) and the (w p i V ot,«i«) plane (right panel), while the bottom row show the 
reconstructed w(z) for the current data (left panel) and future data (right panel). 
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reconstruction applied to inflation. In that prior, about 
92% of models are driven to the r = late-time "attrac- 
tor" solution (using a uniform prior in e) [56]. However, 
in the posterior, the models with tensor/scalar ratio r up 
to the upper limit allowed by the data, which is signif- 
icantly larger than zero, are explored very well [61, 62]. 
In fact, a class of models which has r ~ and a very 
blue spectral index, which forms a fraction of about 90% 
of models generated by the prior, is completely ruled out 
by the data. In this sense, as long as the data have con- 
straining power, the Metropolis-Hastings sampler does 
not necessarily reproduce biases due to the measure on 
the prior which arc seen in a naive Monte-Carlo process. 

Notice however that our thawing/freezing classification 
was initiated at the time when we set initial conditions for 
the scalar field dynamics, z — z s tartj an d did not allow 
the field to settle past the transient and into the late- 
time behavior. If, instead, we require dw/dhia to have 
a given sign only at z < z, where < z s tart, then the 
effect of the initial conditions will not be as important. 
With 2* = 1, and looking at the posterior distribution of 
models, we find that the thawing fraction of models rises 
to 2% - still small, but much bigger than when = z s tart- 

Several examples of dark energy models from our 
chains that are relevant to the thawing/freezing discus- 
sion are shown in Fig. 11. The top panel shows the po- 
tential V(z)/Vq, the middle panel shows the w(z) his- 
tories, and the bottom panel shows the trajectories of 
these models in the w — dw/d\na plane, together with 
shaded regions that should be occupied by thawing and 
freezing models according to Ref. [73]. Clearly, the lat- 
ter plot shows that our models often do not lie solely in 
either region, nor do they necessarily retain their purely 
thawing or freezing behavior. This remains true even if 
we only restrict to very late-time evolution [z < z* = 1, 
say). Fig. 11 also shows a thick (pink) line, a more gen- 
eral lower bound [90] which is obeyed by all monotonic 
quintessence potentials. This latter bound, unlike the 
one from [73], is obeyed by all our models except for 
those where the field samples a section of the potential 
that is not monotonic. Fig. 12 shows the evolution on 
the w - dw/dhia plane for a number of models, with the 
left panel showing the entire evolution from 2 s tart = 3 
to z = 0, and the right panel shows only the evolution 
for z < z* = 1 for the same models. As the various 
examples in these figures illustrate, the disagreement be- 
tween our models and the phenomenologically expected 
bounds from [73] is due to the interplay between the ini- 
tial conditions of the scalar field, the effects of Hubble 
friction, and the shape of the effective potential (specifi- 
cally, non-monotonicity in the second derivative) that our 
more general class of models possesses. For example, the 
short dashed example in Fig. 11 shows freezing behavior 
early on due to the Hubble friction; then the field ve- 
locity increases and a thawing period ensues when H(z) 
has fallen sufficiently; but at the last, most recent stage 
the potential happens to flatten off and the field enters 
the freezing regime again. Note that the more general 
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FIG. 11: Several dark energy histories from our chains that 
are relevant to the thawing/freezing discussion. Top panel: 
the potentials, V(z)/Vo. Middle panel: w(z) histories of 
these models. Bottom panel: trajectories in the w — dw/d In a 
plane, together with shaded regions that should be occupied 
by thawing and freezing models according to Ref. [73] , the or- 
ange diamonds again denoting z — 0. The bottom panel plot 
shows that our models often do not lie solely in either region, 
nor do they necessarily retain their purely thawing or freezing 
behavior. Reasons for the disagreement are discussed in the 
text. The thick (pink) line in this panel shows a more general 
lower bound on monotonic quintessence potentials from [90] 
which is obeyed by all our models satisfying this assumption. 
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FIG. 12: Typical trajectories in the w — dw/dlna plane for a number of models in the posterior, together with shaded regions 
that should be occupied by thawing and freezing models according to Ref. [73], the orange diamonds again denoting 2 = 0. The 
left panel shows the trajectories starting at z s tart = 3 and the right panel shows just the region between 2 = 1 and z = for 
the same models. The thick (pink) line in both panels shows a more general lower bound on monotonic quintessence potentials 
from [90] which is obeyed by all our models satisfying this assumption. 



bound from [90] is again obeyed as long as the sampled 
potential is monotonic. 

Clearly, the thawing/freezing/neither fractions depend 
on the particular priors - the assumed form of V(<p) and 
the initial conditions. Therefore it seems that, at this 
stage when we do not have a good theoretical under- 
standing of dark energy, it is not useful to specialize into 
studying models that are thawing or freezing, because 
their presence or lack thereof is almost certainly not ro- 
bust to fundamental assumptions. Conversely, if we ever 
obtain a precise measurement with a decidedly positive 
or negative value of dw/dlna (or w a , or a%), the thaw- 
ing/freezing picture can help obtain a fuller theoretical 
understanding — a point also made by others advocating 
this approach [73, 89-91]. 



C. How many dark energy parameters? 

More parameters describing the dark energy sector will 
lead to an improved fit to the data, but are they required 
by the data? This question can be answered using the 
Bayesian Information Criterion (BIC) statistic, 

BIC = -21n£ max + £»ln^ data , (21) 

where £ max is the maximum likelihood in a given class of 
models, D is the number of parameters in a model, and 
^Vdata the number of data points used. A BIC difference 
of 2 or more indicates good evidence, and 6 and more 
strong evidence in favor of a model with the smaller BIC. 

Here we compare the ACDM models which have D = 3 
parameters (Ofjjif*, 6 a and il^h 2 ) and the scalar field 
models with two or more additional parameters (so D = 6 



or more) describing the shape of the effective potential 
(ui start , e start and 7y start , plus any higher derivatives of the 
potential). The ACDM model has a slightly higher value 
of the best fit x 2 = — 21n£ max , by about unity (about 
108, compared to about 107 for the two-extra-parameter 
scalar field models, with 119 degrees of freedom). How- 
ever the penalty from the second term in Eq. (21) for the 
scalar field models is large, 3 In 119 « 14.3, and com- 
pletely overwhelms the gain from the improved likeli- 
hood fit. Therefore, the overall BIC evidence (of +13.3) 
is clear: using the current data, there is no compelling 
evidence whatsoever for models more complicated than 
ACDM . A more comprehensive recent analysis using 
Bayesian model selection, albeit using the empirical vari- 
ables (wo, w a ) to describe the evolution of the dark en- 
ergy, can be found in Ref. [93]. 

While ACDM is clearly an excellent fit to the current 
data, there are two important caveats to keep in mind. 
First, it is possible that models different from ones we 
considered have a lower BIC than the ~ 6-parameter 
scalar field models we considered (although it is still 
highly unlikely they will be favored by the BIC criterion) . 
Second, we have shown that the future data will provide 
a much sharper test of dark energy histories, and may 
possibly provide the traction to see a distinct deviation 
from ACDM . Such a development would surely provide 
a much needed hint about the physics of dark energy. 



D. Physical observables 

Since we have obtained constraints using a variety of 
cosmological probes, it is interesting to ask which physi- 
cal quantities (or alternatively, which dynamical aspects 
of our scalar field models) are well constrained by the 
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data. This is a fairly complex question that we largely 
leave for future work — however, before concluding we 
make an interesting observation. We found that distances 
out to intermediate redshifts, z = 1-3, are determined to 
2-3% for the current data. For the future data, the dis- 
tances will be determined to about 0.5%. Therefore, any 
scalar field model that does not preserve the distance to 
rcdshift z <~ 2 to these accuracies will be ruled out. As 
with all other parameter constraints we presented, we can 
expect that these accuracies to be somewhat degraded 
once we allow a completely general dark energy history, 
not just the class of scalar field models with a smooth po- 
tential. Nevertheless, the aforementioned numbers sug- 
gest, for example, that a cosmological probe that has 
ambition to significantly improve upon the current con- 
straints should determine the mean distance to a (single) 
redshift z > 1 to ~ 1% or better. 



VIII. CONCLUSIONS 

The main motivation behind this work was to study 
all possible dark energy histories within a broad class 
of models — chosen here to be the scalar field models 
with an effective potential described by a polynomial se- 
ries. Adapting the Monte Carlo reconstruction formalism 
(previously applied to inflation) to scan a wide range of 
dark energy models, we have generated millions of dark 
energy models and constrained them with the current 
compilation of data from the Supernova Legacy Survey 
[5] , baryon oscillation results from the Sloan Digital Sky 
Survey [71], cosmic microwave background constraints 
from the WMAP experiment [26], along with the mea- 
surement of the Hubble constant from the Hubble Key 
Project [72]. We have also simulated expected future 
data from the same cosmological probes. Instead of at- 
tempting to study comparative merits of specific cosmo- 
logical probes or optimal survey designs, we have concen- 
trated on addressing a more general set of issues regard- 
ing the future quest for dark energy. 

In particular, we have addressed how the theoretical 
prior — working within a specific class of scalar field 
models parameterized with a polynomial series in V(<f>) 
- combines with the cosmological data to produce the 
constraints on cosmological parameters. We found that 
the theoretical prior is significant, as only specific smooth 
dark energy histories are generated, and that their qual- 
itative nature does not vary significantly with the order 
of the polynomial V(<p); see Figs. 1 and 8. Moreover, the 
same theoretical prior excludes certain regions of param- 
eter space (especially in the equation of state parame- 
ters), and makes the posterior distribution highly non- 
Gaussian. As a consequence, the commonly used Fisher 
matrix formalism applied to the equation of state param- 
eters, which assumes Gaussianity of the likelihood, could 
be a terrible approximation to the exact likelihood in a 
specific class of dark energy models, and would signif- 
icantly bias the constraints. A more general approach 



such as the one we have adopted, Markov Chain Monte 
Carlo combined with an exact integration of the dark en- 
ergy evolution equations, is needed. The constraints we 
have obtained on the fundamental parameters are shown 
in Figs. 3 and 4. 

We have computed the principal components of the 
equation of state for our data compilation. As noted 
in Sec. VA, and for reasons listed there, these principal 
components have been pre-computed for the same survey 
but assuming a general DE history; therefore they are 
not uncorrelated. Nevertheless, we find them very use- 
ful — for example, Figure 5 shows that the current data 
roughly imply ct\ < —0.85, thereby imposing an effective 
upper limit to the weighted average of the equation of 
state. Moreover, in Sec. VB we propose a shortcut to 
compute the parameters wq and w a that approximately 
describe the equation of state history w(z) at low rcdshift 
where the constraining power of the data is the greatest. 
Instead of fitting (wa,w a ), we obtain them directly from 
the first two principal components. The left panel of 
Fig. 7 illustrates that the approximate (wo,w a ) descrip- 
tion of the models fits the dark energy histories almost 
perfectly at z < 1 ; we found that even distances to z = 3 
are recovered to within 2-3%. These results suggest that 
converting from the principal components to an arbitrary 
parametrization can be performed with very small biases, 
and essentially at no additional computational cost. 

The constraints in the (wo,w a ) or (w p i V ot,w a ) planes 
indicate that the constraints with future data will have 
the (inverse area) figure of merit improvement of about 
an order of magnitude over the current data. Given 
that our models approach the currently favored ACDM 
cosmology arbitrarily closely, this means that efforts to 
pursue upcoming experimental efforts, such as Planck, 
JDEM, LST, the planned probes of baryon acoustic os- 
cillations, and the improved Hubble constant measure- 
ments are extremely valuable and likely to lead to signif- 
icant improvement in sweeping away a large fraction of 
currently viable models. 

We have also considered the classification of models 
into "thawing" and "freezing" depending on whether 
they are asymptotically receding from or approaching the 
state of zero kinetic energy where the equation of state 
is —1. We found that the thawing and freezing limits 
on monomial-potential scalar field models, presented in 
Ref. [73], are largely not obeyed by our models. As il- 
lustrated in Fig. 11, this is due to the interplay between 
the initial conditions of the scalar field, the effects of 
Hubble friction, and the shape of the effective potential 
(non-monotonicity in the second derivative) that empir- 
ically generated models possess. As we do not have a 
good understanding behind the physical mechanism that 
powers dark energy, we think that the division of models 
into thawing and freezing, or any similar nomenclature, 
is useful only in the context of future experimental con- 
straints: if the constrained phase space ends up favoring 
one of these subclasses, this division can help obtain a 
fuller theoretical understanding 
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Finally note that the approach that we outlined can be 
trivially generalized to consider, for example, the density 
perturbations of dark energy. By solving the perturba- 
tion equations for each model along with the background 
evolution equations as currently done, we can generate 
trajectories using the MCMC approach. Therefore, a ex- 
haustive "scan" through the possible perturbation sce- 
narios can be obtained, paired with the background evo- 
lution histories, and then constrained with cosmological 
data. This approach would nicely complement the previ- 
ous analyses that largely relied on standard phenomeno- 
logical descriptions of the dark energy sector [94-97]. 

This is an excellent time to perform numerically rigor- 
ous, comprehensive analyses of dark energy models be- 
cause the data is finally allowing precision tests of dark 
energy. Not only are the statistical errors respectable 
(see Table II), but most of the experiments are now re- 
porting careful analyses of systematic errors. As we have 
shown in this paper, data expected in the next ~ 10 
years will lead to precision measurements of the param- 
eters describing the effects of dark energy on the expan- 
sion history of the universe. One can sincerely hope, per- 
haps even expect, that these efforts will lead to important 
hints as to the nature and the origin of dark energy. 



models. We use a combination of SNe la, baryon oscilla- 
tion and Hubble constant measurements, together with 
WMAP constraints on the angular size of the sound hori- 
zon, and the physical matter and baryon densities. 

The SNa la data we use are Supernova Legacy Survey 
data from Ref. [5] with N = 115 SNa measurements, 
which include the low-redshift Calan-Tololo sample. The 
observed apparent magnitude is given by 

m, = 5log H d L ( Zi ) + M + a( Si - 1) + (3 Ci + e, , (Al) 

where i = 1 ... 115 runs through the observed SNe. The 
measured values of the apparent magnitude of supernova 
rrii, its rcdshift Zi, stretch Sj, color Ci and total error e» 
are given in the SNLS paper [5] . The nuisance parameter 
Al is a combination of Hubble parameter and absolute 
magnitude of SNe [1]. The likelihood of measuring the 
cosmological parameter set 6 is given by 



£ S Ne(0) = J cM- X 2 /2)T(a)V(p)dMdadp(A2) 
X\0) = *» 2 , (A3) 
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APPENDIX A: CURRENT COSMOLOGICAL 
DATA AND THE LIKELIHOOD FUNCTION 

Here we report on the current cosmological observa- 
tions that we use in order to constrain the cosmological 



where Via) and V(f3) are priors given to the stretch and 
color coefficients. Since M is given a flat prior, the M 
integral can be done analytically [32] to obtain 
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Finally, we follow the SNLS analysis and adopt the 
ACDM values for the stretch and color parameters, a — 
1.52 ± 0.14 and j3 = 1.57 ± 0.15, and marginalize with 
flat priors over the ±l-er values in a and /?. After ex- 
perimenting with several alternative prior widths for the 
stretch and color coefficients, we conclude that their de- 
tails introduce unobservable changes in our results. 

SNe la provide the best constraints at low redshifts 
(z < 1.7). However, Figure 1 shows that there exist 
classes of models that look roughly like the cosmologi- 
cal constant at those redshifts but dominate the energy 
density at z > 2. Those models are not ruled out by the 
SN data, but nevertheless are clearly inconsistent with 
the observed universe as they prevent sufficient structure 
formation. In order to protect against such models, we 
use the information provided by baryon oscillations and 
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CMB measurements of the angular size of the acoustic 
horizon and the physical matter density. More precisely, 
we use the quantity probed by baryon oscillations 



D v (z = 0.35) = 



r{zf 



H{z) 



1/3 



= (1370 ±64) Mpc 
(A7) 

where r(z) is the comoving angular diameter distance to 
redshift z and the measurement comes from the Sloan 
Digital Sky Survey [71]. 

To represent the CMB measurements given by the 
WMAP experiment [25, 26], we do not use the full likeli- 
hood function as its evaluation would be very time con- 
suming for the millions of models we consider. Instead, 
an excellent approximation is to use the single derived 
quantity that is sensitive to dark energy, the angular size 
of the sound horizon [98, 99], which is equal to the ratio 
of the physical size of the sound horizon and the distance 
to recombination 



0a = 



sh 



r(z = 1089) 
where the sound horizon sh is given by 



= (0.595 ±0.002) deg, (A8) 



sh 



V3ff ( 



Jo 



a + Or) 



-1/2 



da, 



(A9) 

where f2 7 and flu are the energy densities in photons 
and radiation (photons plus neutrinos) respectively. The 
quantities f2 7 and flu are given in terms of the CMB 
temperature and are exceedingly accurately measured, 
so we hold them fixed. 

Finally, since 6 a, Q m h 2 and Qbti 2 are our original pa- 
rameters, and we compute the equation of state history 
w(z) for each model, the Hubble constant Hq is neces- 
sarily a derived parameter. To compute the distance to 
decoupling r(z — 1089) we use not the exact w(z), but 
rather the constant effective equation of state w c h, de- 
fined as 



Wcff 



/" 

Jl/(1+1089) 



Ode(«) w(a) da 



(AlO) 



/(1+1089) 



^de(o) da 



The reason we use w e s is that the original WMAP analy- 
sis [26] uses a constant equation of state to arrive at their 
constraints. 

We adopt the the physical baryon density Vt^h 2 as 
a free parameter and give it the prior consistent with 
WMAP, fl b h 2 = 0.0223 ± 0.00074. The physical mat- 
ter density, il m h 2 , is also a free parameter accurately 
constrained by the CMB. Since our models allow more 



complex behavior in the dark energy sector, we conser- 
vatively double the reported ACDM error in this param- 
eter, adopting Q m h 2 = 0.127 ± 0.02. 

Finally we adopt the Hubble Key project measurement 
of the Hubble constant, H = 0.72 ± 0.08 [72]. The final 
likelihood is 



£{0) = £ SNc (0) x £ BAO (0) x £ CMB {0) x £k o (0) (All) 

where the SNe likelihood is specified in Eq. (A3), the 
BAO and Hq are Gaussian with means and standard de- 
viations specified above, and the CMB likelihood consists 
of 0a, Qbh 2 and £l m h 2 constraints, all of which are Gaus- 
sian with means and standard deviations specified above. 



APPENDIX B: FUTURE DATA 

To simulate the future cosmological data, we consider 
the same cosmological probes as for the current data: 
type la supernovae, baryon acoustic oscillations, CMB 
measurements of the sound horizon, Q m h 2 and ttbh 2 , 
and the measurements of the Hubble constant. Since 
we obviously do not know which cosmological model will 
be favored by future data, we center all observables on 
the ACDM model with Q m = 1 - n DE (z = 0) = 0.25, 
w(z) = — 1 and H = 0.72. As described below, we add 
best-guess systematic errors for all future measurements. 

For the SNa data we assume a SNAP-type experiment 
with 2800 SNe distributed in redshift out to z = 1.7 as 
given by the middle curve of Fig. 9 in Ref. [100], and com- 
bined with 300 local supernovae uniformly distributed in 
the z = 0.03 — 0.08 range. We add systematic errors in 
quadrature with intrinsic random Gaussian errors of 0.15 
mag per SNe. The systematic errors create an effective 
error floor of 0.02 (1 + z^)/2.7 mag per bin of Az = 0.1 
centered at redshift Zi. 

Specifications of proposed future baryon oscillation 
surveys and their expected statistical and systematic er- 
rors are described in the Dark Energy Task Force re- 
port [88] (this is partly based on other studies, such as 
[101]). To a very good approximation, baryon oscillations 
provide measurements of the angular diameter distance 
cU(z) and the expansion rate H(z) out to the redshift(s) 
where the source population of galaxies is located. We 
adopt the DETF findings and assume a JDEM-type sur- 
vey covering the redshift range z € [0.5,2.0] of 10,000 
square degrees. For simplicity we assume only three red- 
shift bins centered at z = 0.75, 1.25 and 1.75 with widths 
Az = 0.5 each. With these parameters, formulae in the 
baryon oscillation section of the DETF report indicate 
that the total expected error (statistical plus systematic) 
for a spectroscopic survey is about 1.5% for the measure- 
ments of the angular diameter distance in each bin, and 
2% for measurements of the expansion rate H (z). There- 
fore we assume a total of six BAO observables - three 
distances and three expansion rates - with the aforemen- 
tioned errors. 



21 



Note that, in the DETF error budget, a conservative 
1% systematic error is assumed in each bin for each BAO 
observable. Controlling the systematics to a level lower 
than this, as well as increasing the observed sky area, 
would enable a future JDEM BAO mission to achieve sig- 
nificantly higher sensitivities than assumed here (C. Ben- 
nett, private communication). Similarly, if the JDEM 
SNa systematics can be understood to better than the 
systematics floor we assumed, the sensitivity of that ex- 
periment could be higher. 

For the future CMB experiment we assume Planck 
[102] which is scheduled for launch in 2008. Planck's fidu- 
cial error in £l m h 2 , for a flat ACDM universe, is about 
1%. Since we are considering models with a complicated 



dark energy sector, we conservatively double this uncer- 
tainty, just as in the case with current data. With the 
assumed cosmology, this amounts to assuming Q m h 2 — 
0.1296±0. 00270. The angular size of the acoustic horizon 
for the fiducial cosmology and its expected Planck error 
are 6 A = (0.595 ± 0.00017) deg. Similarly, the baryon 
fraction is assumed to be fl b h 2 = 0.023 ± 0.00015 

Finally, we conservatively assume independent future 
measurements of the Hubble constant will reach the 5% 
level (more accurate measurements are feasible in princi- 
ple, but the systematics remain the primary obstacles). 
Therefore we take H = 0.72 ± 0.036. 

The full likelihood function is given again by Eq. (All). 
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